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Abstract —In this paper, by exploiting the special features 
of temporal correlations of dynamic sparse channels that path 
delays change slowly over time but path gains evolve faster, we 
propose the strnctnred matching pursuit (SMP) algorithm to re¬ 
alize the reconstruction of dynamic sparse channels. Specifically, 
the SMP algorithm divides the path delays of dynamic sparse 
channels into two different parts to be considered separately, i.e., 
the common channel taps and the dynamic channel taps. Based 
on this separation, the proposed SMP algorithm simultaneonsly 
detects the common channel taps of dynamic sparse channels in 
all time slots at first, and then tracks the dynamic channel taps 
in each single time slot individually. Theoretical analysis of the 
proposed SMP algorithm provides a guarantee that the common 
channel taps can be successfully detected with a high probability, 
and the reconstruction distortion of dynamic sparse channels is 
linearly upper bounded by the noise power. Simnlation resnits 
demonstrate that the proposed SMP algorithm has excellent 
reconstruction performance with competitive computational com¬ 
plexity compared with conventional reconstruction algorithms. 

I. Introduction 

In broadband wireless communication systems, channel 
state information is required at the receiver for coherent signal 
detection since wireless channel distorts the received signals, 
especially when the wireless channel is dynamically changing. 
Therefore, accurate channel estimation becomes a fundamental 
problem for broadband wireless communication systems over 
dynamic wireless channels ifTl . 

Various linear channel estimation methods have been pro¬ 
posed in the literature l2l, but their performance can hardly 
meet the simultaneously booming demand of high-rate and 
high-mobility wireless communications. Recently, a lot of 
physical channel measurements have verified that wireless 
channels exhibit sparsity, i.e., the dimension of a wireless 
channel may be large, but the number of channel taps with sig¬ 
nificant power is usually small, especially in broadband wire¬ 
less communication systems a. By exploiting this channel 
sparsity, compressive sensing (CS) algorithms have been used 
to improve the channel reconstruction performance, such as 
orthogonal matching pursuit (OMP) Q, compressive sampling 
matching pursuit (CoSaMP) a, subspace pursuit (SP) a, 
etc. Compared with conventional linear methods, CS-based 
channel reconstruction methods are able to achieve reliable 
channel estimation with reduced training resources 171 . 


Eurther studies have uncovered additional channel charac¬ 
teristics, e.g., the temporal correlations of practical wireless 
channels; path delays change slowly over time, while path 
gains evolve faster lH). By taking these temporal channel 
correlations into account, we have proposed several simultane¬ 
ous CS-based methods to realize simultaneous reconstruction 
of sparse channels, such as adaptive simultaneous OMP (A- 
SOMP) 191 and structured SP (SSP) ITOl algorithms. Different 
from the CS-based methods ignoring the temporal correlations 
of dynamic sparse channels but reconstructing the channel 
in each time slot independently, the simultaneous CS-based 
methods assume that dynamic channels in several consecutive 
time slots share the same path delay set to improve the channel 
estimation performance. However, the path delays of dynamic 
sparse channels may change over time, or even there maybe 
some mutations of the path delays, so the assumption of the 
simultaneous CS-based methods is not always true in practical 
broadband wireless communication systems ii. 

In this paper, we propose a structured CS algorithm called 
structured matching pursuit (SMP) for the reconstruction of 
dynamic sparse channels in broadband wireless communica¬ 
tion systems. By exploiting the temporal correlations of a 
dynamic sparse channel 0, we divide the path delays of the 
dynamic sparse channel into two different parts to be consid¬ 
ered separately, i.e., the common channel taps and the dynamic 
channel taps. Based on this separation, the proposed SMP 
algorithm simultaneously detects the common channel taps of 
the dynamic sparse channel in all time slots at first. Then the 
path delay set of the common channel taps is used as the 
initialization set in the dynamic channel taps tracking process, 
which aims to track the dynamic channel taps and remove 
the fake taps in the initialization set in each time slot. The 
theoretical analysis of the proposed SMP algorithm based on 
the tool of restricted isometry property (RIP) indicates that the 
common channel taps can be successfully detected with a high 
probability, and the reconstruction distortion of the dynamic 
sparse channel is linearly upper bounded by the noise power. 
Numerical simulations show that the proposed SMP algorithm 
with competitive computational complexity has better channel 
reconstruction performance than conventional reconstruction 
algorithms. 

The rest of this paper is organized as follows. Section II 



Fig. 1. Illustration of the dynamic Vehiculai' B channel with a velocity of 
100 km/h in four consecutive time slots. 

presents the system model. Section III addresses the proposed 
SMP algorithm. Section IV presents the performance analysis, 
and simulation results are provided in Section V. Finally, 
conclusions are drawn in Section VI. 

Notation: We use upper-case and lower-case boldface letters 
to denote matrices and vectors, respectively; (•)^, 

(•)^, and II • lip denote the transpose, conjugate transpose, 
matrix inversion, Moore-Penrose matrix inversion, and Ip norm 
operation, respectively; |r| denotes the number of elements in 
set r while ||x||o denotes the number of non-zero elements in 
vector x; hr denotes the entries of the vector h in the set F; 
#r denotes the submatrix comprising the F columns of #. 

II. System Model 

In this paper, we consider the reconstruction of dynamic 
sparse channels H = ‘ ‘ from under¬ 
sampled measurements Y = • • • ,y*-^^] in r time 

slots. The sparsity level, i.e., the maximum number of non¬ 
zero channel taps in is K. The measurement vector 
y^*^ = [yf\y 2 T'’' usually 

obtained through a linear measurement process 

y(‘) = f = (1) 

where ‘ ‘ ‘ denotes the dynamic 

sparse channel in the f-th time slot with N > M, denotes 
the additive white Gaussian noise (AWGN) vector subject to 
the distribution CA/^(0, and r denotes the coherence 

time of the dynamic sparse channel. The measurement matrix 
in ([T]i is usually assumed to be time-invariant to simplify 
the design of communication systems M- 

= $ = An]. ,r. (2) 

Since $ is usually designed in a random way, without loss of 
generality, we assume (t>^cf)^ = 1 and 4>j ~ 0,Vz j ll). 

The temporal correlations of practical wireless channels 
have been verified through analysis and experiments, even 
when the channels are varying fast 121. Fig. [T] illustrates 
the time-domain impulse response of the dynamic Rayleigh 
fading Vehicular B channel with a velocity of lOOkm/h in four 
consecutive time slots, where each time slot corresponds to one 


OFDM symbol duration of 66.7 us in LTE-A systems imi.it is 
clear that path delays of such dynamic sparse channels change 
slowly, while path gains evolve faster. In order to characterize 
the temporal correlations of such dynamic sparse channels, we 
define a common path delay set as 

ro = {i:zer(*),f = l,2,... ,r}, (3) 

where F*^*) = {i : i £ O, ^ 0} denotes the path delay set 
of and O = {1,2, • • • , TV} denotes the entire set. The size 
of the common path delay set is defined as IFqI = L < K, 
where L denotes the temporal correlation degree of dynamic 
sparse channels. Thus we have 

|rW\r(^)| </v-L, (4) 

which means that the number of dynamic channel taps between 
any two time slots is constraint hy K — L. 

III. The Proposed SMP Algorithm 

In this section, we propose the SMP algorithm to reconstruct 
dynamic sparse channels, whose path delays can be divided 
into two parts, i.e., the common channel taps and the dynamic 
channel taps. The key idea of the proposed SMP algorithm 
is that, common channel taps in all r time slots can be 
simultaneously obtained at first, and then the dynamic channel 
taps in each time slot can be individually detected with low 
complexity to refine the final reconstruction. The pseudocode 
of the proposed SMP algorithm is provided in Algorithm [T] 
which is mainly comprised of the following two parts: 


Algorithm 1 Proposed SMP algorithm 

Input: 

Received signals: Y = 

Measurement matrix: $ = [cf)i, 4>2." ' . ^Ar]i 
Maximum sparsity level: K. 

Output: 

Reconstructed channels: H = • • • , 

1: fo = 0, R = Y. 

2 : for A: = 1 to AT do 
3: Z = 

4: fo = Fq U argmax, Aj = l-^i.il- 

5: R = Y - $t Y. 

^0 To 

6: end for 

7: for t = 1 to r do 

8: Z = 1, = $t y, f^ = fg, r^ = y — 

9: while Z = 1 or ||r '||2 < ||r '“^||2 do 

10: Z = Z -f 1. 

11 : f^ = f^“^ U argmaxj 

12: h' = $t,y. 

13: f^ = {K indices with largest magnitude in h*}. 

14: f = y - 4>p,4>t,y. 

15: end while 

16: f = f'-fh(*) = $ty(*). 

17: end for 

18: return H = • • • , 















1) Common Channel Taps Detection (step l^^step 6): The 
first part of the proposed SMP algorithm aims to simulta¬ 
neously estimate the common path delay set Tq^ which is 
expected to contain most of the true common channel taps. 
Similar to our previously proposed A-SOMP algorithm 13 
and SSP algorithm M, we detect the common channel taps 
by iteratively selecting the element that correlates best with the 
residue signal. However, different from the simultaneous CS- 
based methods like A-SOMP l|3 and SSP ifTOll which directly 
use To to reconstruct the dynamic sparse channel, the proposed 
SMP algorithm utilizes Tq as the initialization for the dynamic 
channel taps tracking. 

2) Dynamic Channel Taps Tracking (step y^^step 17): After 
the estimated common path delay set To has been obtained, the 
second part of the proposed SMP algorithm aims to track the 
dynamic channel taps, which are different from one time slot 
to another. In each time slot (we omit the time slot superscript 
t for notational simplicity), the estimated common path delay 
set is used as the initialization of the dynamic tracking process, 
whereby the path delay set T* of size K is maintained but 
refined in each iteration. Specifically, a channel tap, which is 
considered reliable in one iteration but shown to be wrong in 
another, can be added into or removed from the path delay 
set Benefitting from the common channel taps detection, 
the dynamic channel taps tracking process is able to approach 
the true path delay set T rapidly, and this process will be 
terminated when the residue error signal r* is not reduced, 
i.e., ||r '||2 > ||r'-i|| 2 . Finally, the dynamic sparse channel 
can be reconstructed based on the refined path delay set F. 

The proposed SMP algorithm is essentially different from 
the CS-based algorithms which ignore any temporal cor¬ 
relation of dynamic sparse channels a-ini. Meanwhile, it 
also differs from the simultaneous CS-based algorithms 13, 
El, which assume the dynamic sparse channels in several 
consecutive time slots shares exactly the same path delay set. 

IV. Performance Analysis 

To investigate the performance guarantee of the proposed 
SMP algorithm, in this section, we first derive the condition 
of the successful detection of the common channel taps, and 
then analyze the condition of the successful tracking of the 
dynamic channel taps. 

A widely used condition of $ for reconstruction perfor¬ 
mance analysis is named RIP lb), whereby a sensing matrix 
$ is said to satisfy the RIP of order K if there exists a constant 
5 G (0,1) such that 

(l-b)||x||2<||$x||2<(l + J)||x||^ (5) 

holds true for all AT-sparse vectors x G (||x||o = K). 

In particular, the minimum of all constant 5 satisfying Q is 
referred to as a restricted isometry constant 5k- 

A. Common Channel Taps Detection 

For the common channel taps detection, the proposed SMP 
algorithm selects the element that maximizes the sum of 
absolute correlation between the basis of $ and the residue 


measurement R. This greedy selection criterion can be for¬ 
malized as follows: 

T 

argmaxAi = ^ i = 1, 2, • • • , (6) 

* 

where Zi t denotes the correlation between the i-th basis of $ 
and the residue measurement in the f-th time slot. 

In the first iteration, which means that the residue measure¬ 
ment R equals to Y and the rough common path delay set Fo 
is empty, the index i is selected if and only if 

Xi > Xj, Vj G H and j ^ i. (7) 

Here we introduce two results of the RIP condition (|5]i. 
The first one is named near-orthogonality 161: Let /, J C H 
be two disjoint sets, i.e., / fl J = 0, if b| 7 |+|j| < 1, then for 
an arbitrary vector x G we have 

||$f$,/x||2 < (5 |/|+|J|||x||2. (8) 

The second one is called RIP inequality For an arbitrary 
vector X G where / C H, we have 

(1 - ^|7|)||x|| 2 < ll^f ^7 x||2 < (1 + 5|7|)I|x||2, (9) 

--^||X||2 < ||($f$7)-'x||2 < r^ll^ll2- (10) 

1 -b 0|/| 1 - d|/| 

Then, by applying near-orthogonality and RIP inequality, 
Xi in (|6l) can be derived as 

A, = ^|c/.fy0)|=^|c/.f($h(‘)+n0))| 

> ^(Ic/.f ^ - I 0 f n 0 )|) 

t=l 

> E(l#|-E^2|/ii0|-|c/,fn0)|). (11) 

t—1 s^i 

Similarly, Xj can be derived as 

A, < I + E^2|/iiO| + |c/.fn0)|). (12) 

t=l siCj 

Using (|7]i, (fTTI) . and (fT2l) . we can obtain the sufficient condition 
under which the index i will be selected: 

r r 

E ^ I ^ E I + “U '^j G ^ and j ^ i, (13) 

where the parameter aj is 

= E(E^2|/iiO| +E^2|/iiO| + |c/,fn0)| + |c/.f n0)|). 

s/j s^i 

(14) 

From (O, it is clear that when the sum of the absolute 
values of a channel tap is larger than that of other channel 
taps plus a small value ctj, the index of this channel tap will 
be added into the rough common path delay set Fq- The size 
of the rough common path delay set Fq is K, while there are 




only L < K common channel taps. Thus, the true L common 
channel taps will be contained in Tq with a great probability. 
As shown in Fig. [T] the sum of the absolute values of a 
common channel tap is much larger than that of the dynamic 
channel taps. This is because the absolute values of a common 
channel tap in different time slots can be summed together 
while not true for the dynamic channel taps. 


B. Dynamic Channel Taps Tracking 

Intuitively, in each time slot (we omit the time slot super¬ 
script t for notational simplicity), the proposed SMP algorithm 
tracks the dynamic channel taps to refine the final reconstruc¬ 
tion, whereby a reliable channel tap will be added into and an 
unreliable one will be removed from the estimated path delay 
set r in each iteration. In order to measure the performance of 
the dynamic channel taps tracking process, the reconstruction 
distortion ||h — h ||2 is important to be analyzed. 

Here we introduce the residue-orthogonality @ of the 
projection operator used in the CS theory: For an arbitrary 
vector X € and a sampling matrix S (jArx|/| 
column rank, we have 

J $fresid(x,$/) = 0, 

I l|resid(x ,$/)||2 < ||x|| 2 , 


where resid(x, $/) = x — $/$|x denotes the residue signal. 

Thus, by applying RIP condition and residue-orthogonality, 
||r ^||2 of the proposed SMP algorithm can be derived as 

I|r'll2 = ||resid(y,$p,)ll2 

< ||resid($p^p,hp^p,,$pi)ll 2 + ||resid(n, $p ,)||2 

— Il^r\f‘^r\f'II 2 + 

< \/l + iJ2if+i||hp,^pi II 2 + ||n||2. (16) 

Similarly, ||r *~^||2 can be also derived as 

||r'-i|l2 = ||resid(y,$p,_0ll2 

> l=i^^||hp^p,_,||2-||n||2. (17) 

In order to compare ||r *||2 and ||r^“^|| 2 , we need to obtain 
the relationship between ||hp^pi ||2 and ||hp^pi_i II 2 . From the 
set F* in the process of dynamic channel taps tracking, we can 
derive the following two inequalities 


||hr\pi ||2 ^ 262K+1 

\/K “ (1 — (52if-i-i) 


IIhp^pi II 2 ^ 


1 + 52K+1 

1 — ^2K+1 


2 ll^r\f*-i II 

||hpypi II 2 + 


2(1 -I- 52k+i) I 

1 — ^2K+l 

2 „ „ 


1 — ^2K+1 


n 2 . 


(18) 

(19) 


The proof of these two inequalities is similar to that of 
Appendices H and I in ||6l, and we omit the proof here due 
to their similarity and paper space limitation. By applying 
dUll and (EUl into (QSll and (Ell, the termination condition 
I|r'll 2 > ||r'-i ||2 for dynamic channel taps tracking can lead 
to the result as 


|n ||2 > 


52K+1 I 

ViT + 2 




^2if-|-l I 

/^+2 


hr\fll2, (20) 


which means that if the noise power is larger than that of the 
residue channel taps multiplied with a coefficient, the dynamic 
tracking process will be terminated, and the final estimated 
path delay set can be obtained as F = F^”^. 

Then we analyze the relationship between the reconstruction 
distortion ||h— h ||2 and ||hp^p ||2 as 


|h-h||2 < ||hp - $ty||2 + ||hp^p||2 


< ||hp-$I($rhr)||2 


+ ll^pn||2 + 


^r\f II 


< 


1 — S 2 K+I 1 ~ 52K+1 


I ^ ^2K-\-l II II 1 \ 

2 + -?-l|n|| 2 , (21) 


where hp = €>p^>ty and h^^p = 0. By substituting dTlI i into 
®, we have 


Ih-hIU < 


2 -f ^/K -(- S 2 K +1 + I 

52K+i{^ — ^2K+l) 


( 22 ) 


which means that the reconstruction distortion is linearly upper 
bounded by the noise power ||n|| 2 , the sparsity level K, and 
the RIP of the measurement matrix $. 

In addition, the main computational burden of the proposed 
SMP algorithm comes from the matrix pseudo inversion calcu¬ 
lating, which can be effectively solved by the Gram-Schmidt 
algorithm na. Compared with the conventional reconstruction 
algorithms, the complexity of the proposed SMP algorithm is 
competitive. 


V. Simulation Results 


In this section, we investigate the performance of the 
proposed SMP algorithm through numerical simulations. We 
construct the measurement matrix $ of size M x N = 
100 X 200 whose elements are independent and identically 
distributed (i.i.d.), i.e., (j)ij ^ Af{0,l/M). The gains of 
channel taps at each time slot are drawn from the standard 
Gaussian distribution, i.e., ^ A/’(0, l),i S while the 

path delays of the dynamic channel taps for the f-th time slot 
are randomly chosen from H \ Fq. The coherence time of the 
dynamic sparse channel is set as t = 10. 

A commonly used metric to evaluate the channel recon¬ 
struction performance in practical communication systems is 
the mean squared error (MSE) defined as 


MSE = ^ ^ 


||h(‘) -hW||2 

||h(‘)||2 


(23) 


In the simulation, we perform 10,000 independent trials to 
obtain the MSE performance of the following reconstruction 
algorithms: 1) Linear method ||2l; 2) OMP HD; 3) SP 0; 
4) A-SOMP ||9l; 5) Oracle least square (oracle-LS) ifTSl . 
Specifically, the oracle-LS is regarded as the theoretical bound 
of the MSE performance due to the fact that the path delay 
set is assumed to be perfectly known as a prior information. 

Eig.|2]investigates the MSE performance comparison against 
SNR for different reconstruction algorithms, where the tem¬ 
poral correlation degree is set as L = 5. It is evident that the 
standard OMP algorithm and the SP algorithm outperform the 

















Fig. 2. MSE performance comparison against the SNR. 



linear method by about 2 dB, where the benefit comes from 
utilizing the channel sparsity. Further, the A-SOMP algorithm 
is better than the SP algorithm by about 2 dB when the SNR 
is low, since it roughly considers the temporal correlations of 
the dynamic sparse channel. For the proposed SMP algorithm, 
it is clear that another 2 dB SNR gain can be achieved due 
to its capability to detect the common channel taps accurately 
and track the dynamic channel taps rapidly as discussed in 
Section IV. 

Fig. m shows the MSE performance comparison against 
the temporal correlation degree L, where SNR = 5 dB 
is considered. It is clear that the SP algorithm achieves 
lower error levels than the standard OMP algorithm, while 
the conventional linear method performs worst. The MSE 
performance of the A-SOMP algorithm is competitive only 
when the temporal correlation degree L is very large, while 
the MSE performance of the proposed SMP algorithm is the 
best for all considered temporal correlation degree. It should be 
pointed out that when there are no dynamic channel taps, i.e., 
L = K, the MSE performance of the proposed SMP algorithm 
and the A-SOMP algorithm achieves the performance bound 
of the oracle-LS method. 


VI. Conclusions 

In this paper, we proposed a structured CS algorithm called 
SMP to reconstruct dynamic sparse channels for broadband 
wireless communication systems. By exploiting the temporal 
correlations of dynamic sparse channels, the path delays are 
divided into two parts, i.e., the common channel taps and the 
dynamic channel taps. Based on this separation, the common 
channel taps over all time slots are simultaneously detected at 
first. Then, the proposed SMP algorithm tracks the dynamic 
channel taps for each single time slot individually. The theo¬ 
retical analysis of the proposed SMP algorithm provides the 
condition of successfully detecting the common channel taps 
as well as the upper bound of the reconstruction distortion. 
Einally, simulation results demonstrated that the proposed 
SMP algorithm is able to accurately reconstruct dynamic 
sparse channels with about 2 dB SNR gain compared with 
the recently proposed simultaneous CS-based reconstruction 
algorithms. 
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